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Abstract 

An unstable gravitino with lifetime longer than 10^^ sec or so has been proposed as a possible dark matter 
candidate in supergravity models with R-parity breaking. We find a natural realization of this idea in the 
minimal supersymmetric left-right models where left-right symmetry breaking scale in the few TeV range. 
It is known that in these models, R-parity must break in order to have parity breaking as required by low 
energy weak interactions. The sub-eV neutrino masses imply that R-parity breaking effects in this model 
must be highly suppressed. This in turn makes the gravitino LSP long lived enough, so that it becomes the 
dark matter of the Universe. It also allows detectable displaced vertices at the LHC from NLSP decays. We 
present a detailed analysis of the model and some aspects of its rich phenomenology. 
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I. INTRODUCTION 



It is now widely believed that about 25% of the energy density in the Universe is in the form of 
a cold dark matter. The nature of the elementary particle which constitutes this dark component is 
however not known. The scenario with an unstable gravitino which is the lightest supersymmetric 
particle (LSP) with a very long life time as a dark matter is particularly attractive lli |20, having 
interesting consequences for astrophysics and cosmology, as well as colliders. 

The LSP gravitino in an R-parity conserving supersymmetric (SUSY) theory is absolutely sta- 
ble and could be a dark matter candidate M\. An unstable gravitino, in a theory with broken 
R-parity symmetry, needs to be sufficiently long lived to be a dark matter. The strengths of the 
R-parity violating (RPV) couplings (usually denoted by A, A') needed to achive required longevity 
must be highly suppressed i.e. A, A' < 10^^, and one would then like to understand the origin of 
such small couplings. It is therefore interesting to explore models where such small couplings may 
arise naturally. We find that TeV scale supersymmetric left-right (SUSYLR) models which have 
been discussed in connection with neutrino masses and strong CP problems provides one such 
framework, with two interesting features: (a) R-parity breaking is dynamically induced in the 
global minimum of the potential in order for the theory to break parity; (b) smallness of neutrino 
masses guarantees that the resulting R-parity violating interactions are highly suppressed. 

The left-right symmetric theories based on the gauge group SU (2)^ x SU (2)/j x U (1) b-lM 

were originally introduced to explain parity violation in the standard model and were found to have 

a number of interesting properties. These models provide a natural framework for understanding 

the small neutrino masses via the seesaw mechanism[5] and became specially interesting after the 

discovery of neutrino masses. The right handed neutrino which is essential for seesaw mechanism 

n 

is automatically contained and furthermore the B-L gauge symmetry[6] whose breaking provides 
the heavy Majorana mass to the right handed neutrinos relates the small neutrino masses to the 
parity breaking scale. This implementation of the seesaw mechanism is different from those in the 
SM context where the right handed neutrino mass is a free parameter. 

The seesaw scale (or the parity breaking scale) however, still remains undetermined. It can 
take values anywhere from near a TeV if the Dirac Yukawa couplings ~ 10^® ~ or it can 
be much higher if the Y„'s are larger. It is important to note that a TeV seesaw scale is perfectly 
natural and does not require any higher fine tuning than that present in SM. A TeV seesaw scale is 
clearly of great interest for the LHC. It is worth noting that the current searches by the CDF and 
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DO collaborations have yielded limits on the parity breaking scale in the 750 GeV range ||7|]. 

In this paper we consider the minimal left-right (SUSYLR) extension of the MSSM since it 
combines the advantages of supersymmetry while making MSSM realistic by providing a way 
to understand neutrino masses. It also cures certain problems of the MSSM such as making R- 
parity a good symmetry if one uses B-L=2 Higgs triplets to break parity so that proton is stable. 
Furthermore, it provides a solution to the SUSY and strong CP problems. The minimal version of 
this model (without any gauge singlets) has two striking features: (i) the parity and SUSY breaking 
scales are related, thereby predicting that the seesaw scale is necessarily in the TeV range and (ii) 
while R-parity conservation is automatic above the parity breaking scale, the ground state of the 
theory can break parity only if R-parity is spontaneously broken by a vaccum expectation value 
(VEV) of the right-handed sneutrino fields [8]. Thus the model breaks R-parity dynamically. 
Furthermore in this case there is an upper limit on the Wr scale in the range of a few TeV's 
which is to be expected since parity breaking and susy breaking are intimately linked. This makes 
it possible to test the theory using LHC data expected in coming years. 

An interesting feature of the resulting dynamical R-parity breaking is that it only breaks lepton 
number and keeps baryon number intact and therefore the proton is absolutely stable in this model. 
It is also worth noting that while the effective R-parity breaking below the Wji scale has some 
properties similar to MSSM with bilinear R-parity breaking [|l2n . it has many properties which are 
characteristic of the SUSYLR theory that can provide distinguishing tests. 

The immediate question that then arises is whether this SUSY model which starts out promising 
a stable dark matter, does indeed have a dark matter after parity and R-parity breaking. We address 
this question in this paper. We find that despite R-parity breaking, the unstable LSP gravitino in 
our model can be the dark matter of the universe [1]. The reason for this is that requiring sub- 
eV neutrino masses suppresses the strengths of the R-parity breaking interactions responsible for 
gravitino decay to such a level that the gravitino becomes sufficiently long lived and dark matter. 
Secondly, we also find that despite the neutrino-Higgsino mixing induced by spontaneous R-parity 
breaking, the seesaw results for neutrino masses remain essentially intact. Another consequence 
of this small strength of R-parity breaking is that the next-to-lightest SUSY particle (NLSP) which 
can be a neutralino or stau or sneutrino produced at LHC has a lifetime such that it can give rise 
to displaced vertices in the LHC detector [13]. 

The paper is organized as follows: Section II presents the basic contents of the model, and 
reviews the result of Ref. [|80 that R-parity indeed must break in the model if ordinary parity has 
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to break. In Section in, we consider implications of RPV on the neutrino mass in this model. In 
Section IV, we discuss contribution to various RPV couplings due to bilinear terms in the super- 
potential. In Section V, we derive cosmological implications for gravitino as the dark matter. In 

Section VI, we discuss collider signatures of the model and in Section VII, we briefly discuss some 
other consequences of the model. In the appendix, we display the minimization of the potential 
and obtain bounds on the H^ij-boson mass and the right-handed sneutrino VEV 

n. BASIC FEATURES OF SUSYLR MODEL AND SPONTANEOUS R-PARITY VIOLATION 

The gauge group in the minimal SUSYLR model is SU{2)l x SU{2)r x U{1)b-l x 
SU{3)c. The chiral left and right handed quark superfields are Q = (m, rf)(2, 1, |, 3) and 
Q'^ = (m'^, 2, — |, 3*) respectively, and similarly the lepton superfields are given by L = 
{u, e)(2, 1, —1, 1) and L'^ = (u'^, e'^)(2, 1, +1, 1), where flavor indices have been implicit. The 
symmetry breaking is achieved by the following bi-fundamental and B — L — ±2 triplet Higgs su- 
perfields: 0„(2, 2, 0, 1) (a = 1, 2), A(3, 1, +2, 1), A(3, 1, -2, 1), A^(3, 1, -2, 1), A^(3, 1, +2, 1). 
The superpotential of the model is: 

+ y,L^r2$ir2L^ + y5L^r2$2T-2L^ 

+ if {L^T2^L + L'=^T2A'=L'=) + HabTx {^It^^t^) 

+ //aTt (AA + A^A^) , (1) 

where F's are Yukawa couplings, / is the Majorana coupling and //a is the /x-term for triplets. 
Note that we do not have any gauge singlet fields in the model. 

First point to note is that if Mr » Msusy, the SUSY breaking scale, the right handed gauge 
symmetry remains unbroken since we must have (A*^) = (A*^) = to preserve supersymmetry. 
From this it follows that these vevs i.e. (A'^) = vr and (A*^) = vr and yU-term must have TeV 
scale vevs, i.e., the theory is necessarily one with TeV scale parity breaking. From now on we will 
assume that all the mass parameters in the theory are of order of a TeV. 

To do a detailed analysis of the ground state of the theory, we write down the scalar potential 
including the soft SUSY-breaking terms: 

V ^Vf + Vd+Vs, (2) 
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We first give Vd 

Vd 



9l 



J2 li^'^i^'S^.s + 2Tr (A^t^^^A^ + A'^r^A') + Tr($r^<l>t)| 

m 

2 

+ |^^|2Tr(AV^A + AV„A)+Tr($V„<l>)|' 

m 

+ ^ |i>^t^>c ^ 2Tr {A^A - A^^A' + A+A - A^^^^) _ 



(3) 



Below we give some of the terms from Vp = J^, ^ with 0^ going over all the fields in the 



model and the soft susy breaking term Vs M, 



14|] which are relevant for our discussion: 



Vf + Vs 



(4) 



+ iA (j7t2AL + L'^^TsA'^L" 

+ hab^V {^lT2^bT2) + BTl (AA + AW) + h.c] 

+ M2,Tr($t$,)+m^(LtL + L'=tL^) 

+ MiTr (A"fA + A^^A^) + M|Tr (A+A + A^^^^) ^ 

where in addition to the terms dependent on Higgs fields, we have also kept slepton terms since 
sneutrino is electrically neutral and can in principle have vev. We have omitted all terms involving 
the squarks. All soft parameters with mass dimensions are assumed to have TeV SUSY-breaking 
scale. 

It was shown in [|81] that if the ground state of this tree level potential has to break parity, it must 
break R-parity by giving a vev to the field. We review this argument in the appendix. 

The global minimum of the model is then characterized by the following vev pattern of the 
fields 




In this case, there is an upper bound on the vr scale as discussed 

mo (A + //iA) 




(5) 



given roughly by 



(6) 



where the parameters A and mo are susy breaking parameters in the theory, of the order of a few 
hundred GeVs and / is a typical Majorana coupling of the A fields. Using these, we get few TeV 
upper limit on the right-handed scale if / ~ ^. A detailed derivation of the upper bound is given 
in the appendix. 



It has recently been shown that in the case of the model with an additional singlet, once one 
takes one loop corrections into account, in the domain of parameters where m| < 0, R-parity 
conserving and electric charge conserving minimum indeed becomes the global minimum liol]. In 
this case one generically needs a higher parity breaking scale to get correct mass spectrum for 
sleptons. Our discussions in the subsequent sections will remain valid for the complementary 
domain i.e. m| > in the presence of loop effects. R-parity conserving ground state can also 
be a global minimum if one includes higher dimensional operators provided the parity breaking 
scale is above 10^° GeV llllll . Since we are interested in TeV scale Wr, these new operators do not 
affect our considerations. 

In the rest of the paper, we will assume that the component of the right handed sneutrinos that 
acquires a vev must align along the electron flavor direction so that the neutrino masses will remain 
in the sub-eV range. It appears that it may be possible to ensure this by choosing the Aee term 
associated with u^u^A'^ coupling to be sufficiently negative while keeping this for other flavors to 
be positive. We do not pursue the details of this calculation here. 



III. NEUTRINO MASS 



In this section, we address the question of how to understand small neutrino masses in this 
model. As is well known, neutrinos acquire Dirac masses after electroweak symmetry is broken 
by the $1,2 vevs. We assume that only < > and < > fields acquire vevs ki and K2 re- 
spectively. In this model, B-L breaking gives large Majorana masses to the right handed neutrinos 
which combined with the the Dirac masses leads to the usual type I seesaw mechanism. In the 
absence of the sneutrino vev, the type I seesaw mechanism yields sub-eV left-handed neutrino 
masses from TeV-scale right handed-neutrino masses provided the neutrino Dirac Yukawa cou- 
plings are of the same order as that of the electron in the standard model. An upper bound on 
the vr scale dieted by the dynamics of the model, implies that all elements of the neutrino Yukawa 
coupling matrix in our model must have an upper bound of order 10^. There are no type n 
contributions in the renormalizable supersymmetric left-right model in the supersymmetric limit. 
After SUSY breaking, a finite but small type II contribution is induced which can be comparable 
to the type I contribution. This does not affect the analysis of neutrino masses done below and 
we ignore it here. As we show later, the small Dirac Yukawa couplings needed for understanding 
small neutrino masses have the important implication, that the gravitino lifetime which is inversely 
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proportional to is long enough that it can become the dark matter of the Universe. 

The presence of right handed sneutrino condensate complicates the analysis of neutrino masses 
since it introduces a mixing term of the form LH^ in the superpotential making the neutrino- 
Higgsino-gaugino mass matrix to be a 15 x 15 mass matrix. As we show below, if we assume that 
the right handed sneutrino that acquires a vev is aligned along the electron flavor direction, one 
can still use successive seesaw approximation so that neutrino masses remain in the sub-eV range 
as in the simple type I seesaw models prior to neutrino Higgsino mixing. 

To proceed with the neutrino masses, note that the terms in the superpotential are induced by 
the sneutrino vev and expand the seesaw matrix are given by 

AW = iY,)ijXjLj{iT2)LH^ + {Yi)ijXjLj{iT2)LK , (7) 

where H^'s are defined through $i = {H'^,Hy) and $2 = {Hd^H'^), and Xi — {vf). These 
terms mix lepton superfields with the Higgs superfields and therefore break R-parity. They lead 
to the mixing of lefthanded neutrinos with the neutralinos. These mixing terms enlarge the 6x6 
normal seesaw neutrino mass matrix to 15 x 15. However, we can employ successive decoupling to 
simplify the analysis of this matrix to get neutrino masses. For f ^ in the multi-TeV range, the v^- 
scale fermions decouple first, leaving a 9 x 9 mass matrix. On the further application of a second 
stage seesaw through integrating out the electroweak scale fields, small neutrino masses follow. 
For convenience of discussion, we choose a basis in which the Yi is diagonal and all neutrino 
mixings arise from Y^,. We will also choose all Majorana couplings of right handed neutrinos, /, 
to be diagonal. 

As noted before we will assume the RH-sneutrino vev to align along the elec- 
tron flavor i.e. (^^;E1,t) = so that the neutrino-Higgsino mixing is in the MeV 
range. The symmetric neutrino-neutralino mass matrix can now be written in the basis of 
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aBL),i^f.AO^A"'^} as 
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where we have assumed the vev's of doublet fields have the following pattern 

$1 





0^ 




f K2 
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Ki j 




I 





(8) 



To simplify this matrix to obtain the neutrino masses, we proceed in two steps as mentioned: first 
we integrate out the TeV scale fields uf, A'^, A'^ and a combination of X%, Abl- The remaining 
combination B = , „ (—iX%) H — ^J&=(— zAbl) stays at the electroweak scale. The mass 



BL 



matrix for the electroweak scale active fields is 
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B 



where 
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(9) 



-gi^^R + ^^t^Mbl and ^' satisfies ^'-^ = g-^ + g'l 
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When Yi, and Yi are small, we can still use the seesaw formula to estimate the neutrino mass 
contributions from R-parity violating terms. 



6M^ 



[(XuhkXk, 0, {Yi)ikXk, 0, 0, 0] 
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Ml 










- 9' K-1 













Mb 








(10) 



The sneutrino condensates picks up 11, 13, 31 and 33 elements of the inverse matrix. To apply 
the seesaw formula, we concentrate on the corresponding cofactors of these elements. 



9l 



+ 



M, 



BL 



C0fi3 = COf: 



gla'^ 2 

31 — 7\ — /^12/^ll/t2 

2 V 91 



Ml Mr Mbl 

2 ^ ^2 + 



9r 



91l 



glg'^ 2 2 (Ml 

COfgs = /^11«2 ( 



Mr M, 



BL 



9l 



91l 



(11) 



Therefore, the neutrino mass can be written as 



{MA 



- !^(Y f-^Y^) 

VR^ "-^ 2MLM^(/.n//22-/^^ 



9l9'^4 



I2) 



Ml Mr Mbl 

2 „2 



9l 9r 



9bl 



X 



— ) {yu)ik{yu)jk + xl{Y^)ayeSji + xl{Y^)ijyeSii - ( — ) y^i^ji 

1^11/ \/^12/ 



(12) 



where the first term is the usual see-saw formula for neutrino; we have worked in the basis that the 
charged lepton Yukawa coupling is diagonal and assumed only xi ^ 0. In the following sections, 
we will see that this constraint helps us to avoid too large RPV terms, which not only makes 
gravitino life-time sufficiently long but also has important collider implications. 

An interesting aspect of the neutrino mass formula is the presence of the term 
In the framework of mSUGRA and gauge mediated SUSY breaking, one 



91 ~^ 9% 



"sir 



typically gets 



\Ml 
9l 



IM, 



R\ 



IM, 



BL| 



9\ 



9L 



(13) 



This relation is preserved by renormalization group evolution. In SUSYLR case, we have Ml — 
M^ can be complex, while Mbl must be real. In a special case where arg(ML) = ±27r/3, the 



second term in Eq. (fT2l) vanishes, i.e. 5 My = 0. In this case, we can also allow X2, 7^ without 
affecting the neutrino mass discussion above. However no such assumption is necessary if we only 
choose Xi 7^ 0. 

Note that in the context of MSSM with R-parity breaking bilinear terms, getting small neu- 



23D. In 



trino masses generically requires to fine-tune the neutrino-Higgsino coupling terms 1221 
our case, the smallness of these terms is now guaranteed by the smallness of the Dirac Yukawa 
couplings together with the alignment of {u'^) along the electron direction. 



IV. SOURCES OF R-PARITY VIOLATION IN THE MODEL 



Our model is different from all previous models of R-parity breaking [|2 11 1221] because R-parity 
breaking is forced by the tree level dynamics of the theory to enable parity breaking [8]: We do 
not have the freedom of choosing the tree level parameters leading to an alternative vacuum with 
R-parity conservation. Here R-parity is broken by the right-handed sneutrino vev which leads to 
bilinear R-parity violating coupling at low energies. As noted recently, in models with singlets, 
once one loop corrections are included lioi . in a parameter subdomain, one can indeed have an 
R-parity conserving vacuum as a global minimum without breaking electric charge. If the same 
discussion were to apply to the minimal case under discussion here, the considerations of this 
paper will still remain valid in the complementary parameter domain. 

The dominant R-parity violating interaction in the effective theory below the TeV scale, after 
the right-handed sneutrino and the heavy A fields are integrated out, is given by 

AW = e,Lf{iT2)H^ + e[mtT2)K ■ (14) 
There are analogous terms coming from the soft-breaking trilinear terms given by 

Ci = B,Lj{tT2)H^ + B[Lj{tT2)K, (15) 

where Cj = {Yy)iiXi and e[ = YeXiSn. The terms in the second equation above arise when we 
implement the soft SUSY breaking on the effective R-parity breaking Lagrangian (say via usual 
mSUGRA). Since in general the coefficients are arbitrary and we ignore then in what follows. 

In order to understand neutrino masses, we have assumed that only the right-handed electron 
sneutrino gets a VEV. The main consequence of this is that since Yg ~ 10^^^ and Y^ ij < 10~^, 
the strengths of RPV interactions are all very small. In the presence of these terms, the gravitino is 
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unstable; however, we will see later that the small magnitude of these coupling strengths required 
to understand neutrino masses, allows the gravitino to live long enough to be the dark matter of 
the Universe. Thus the gravitino dark matter is connected to the smallness of the neutrino mass, 
and the viability of the scenario does not require tuning of a separate parameter. 

I 




Z (W) 

FIG. 1: "Flavor- violating" gauge couplings due to bilinear R-parity violation. The shaded blob is an R- 
Parity breaking vertex. 



<v^> 



H 



H 



FIG. 2: The induced tri-Unear RPV couplings. The shaded blobs are R-Parity breaking vertices. 



The bilinear R-parity breaking terms in the superpotential lead to mixings between neutralinos, 
Hu and neutrinos u as discussed in the previous section and similarly between charginos and 
charged leptons. These terms can also generate mixing between sleptons and boson through 
the F-term of Hu and H'^ 

Fhu = /^ii-f^d + + Q-^i , 

Fh'^ = A*2i^^d + li22Hd + e'lLi . (16) 
11 



To calculate R-parity violating decays, a general way is to diagonalize all the mass matrices and 
then to rotate all fields into physical states. In the present work, since we only estimate the order 
of magnitude for the R-parity violating effects, we adopt an approximate but more convenient 
approach. We take the bilinear RPV terms as mass perturbations which are to be inserted into 
all the amplitudes. It is these mass insertions that act as sources of R-parity violation. In Fig^ 



we illustrate how to get the "flavor- violating" gauge coupling induced by bilinear RPV term [12711 . 
which plays an important role in neutralino decay. The tri-linear couplings of the conventional 
R-parity breaking XLLE^ and X'QLD^ terms can also be generated from the bilinear terms. The 
diagram on the left-panel in Fig. 2 is generated from neutrino-neutralino mixing in Eq. dH), and 
that on the right is generated through slepton-Higgs mixing due to Eq. (fT6l ). As expected, the 
right-handed (s)leptons can be coupled to leptons, whereas the left-handed (s)leptons couple to 
both hadrons and leptons. 

V, GRAVITINO DARK MATTER 

In generic SUGRA theories, the gravitino is a very weakly coupled particle with mass ranging 
from eV to many TeV's. It can be produced in early universe plasma and remains in equilibrium 
with rest of the cosmic soup at very high twemperatures i.e. (T ~ Mpi). Slightly once the universe 
cools below the Planck temperature, the gravitinos decouple. Since their annihilation or decay rate 
are very slow, their number density dilutes only due to entropy dumped into the cosmic bath at 
different annihilation thresholds of other particles. This dilution is not a large effect. Therefore 
in the absence of inflation, if gravitino is the LSP and R-parity is conserved, its mass must not to 
exceed 1 keV in order not to over-close the universe. 

In the inflationary scenario however, any initial gravitino abundance will be diluted to very tiny 
values. However, secondary production of gravitinos in the reheating process can be appreciable 
and proportional to the reheat temperature. This has been estimated in various papers to be jisl . 



where Tr is the reheating temperature, and m3/2 and rrig are the gravitino and gluino masses, re- 
spectively. While this would permit gravitinos in the 100 GeV mass range, one must be mindful 
of other constraints when R-parity is conserved: the NLSP (often neutralino or stau, etc.) decays 
late, after freeze-out and often after big-bang nucleosynthesis (BBN) (T < 1 MeV), and produces 
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a large amount of entropy mainly as photons into the universe, thereby drastically changing the 
ratio of — . Also the decay products can destroy the produced elements making it hard to under- 
stand the successes of BBN [16]. Therefore, in a consistent picture of the universe described by 
supergravity theories, the NLSP must decay quickly (< lO^s). If R-pairty is conserved, typical 
NLSP life time however is anywhere from a few days to years as can be inferred from the formula 

/ m?j2 \2 /100GeV\^ 



''"NLSP 



V ^NLSP / 

Here we need ^ 1 GeV for mNLSP mass around 100 GeV. So whereas the overclosure 

constraints by gravitinos can be reconciled with multi-GeV mass gravitinos by adjusting the reheat 
temperature, the NLSP lifetime constraints cannot be accommodated without extreme fine tuning. 

This late-decay problem can be solved if R-parity is violated by a small amount yj]. [For an 
R-parity conserving theory, there are also ways to avoid this; see jnll .l In this case, the neutralino 
decays very quickly into RPV channels and there is no longer an upper-bound from BBN consid- 
erations. Our model leads naturally to this scenario since the gravitino can now decay to several 
channels, as shown in Fig. 3. To see if it lives long enough to become a viable dark matter, we 
need to estimate its lifetime. 




G "MXPJXj 



<v^> 




q 



<v^> <^ 



FIG. 3: Left-panel: Dominant PRV two-body gravitino decay mode {G in our model. Right-panel: 

a typical three-body gravitino decay mode. 



The dominant decay-modes of gravitino are 7 + z/, + + u and + z/ [ |18l ]. From the 

following analysis of cosmological constraints, we find the gravitino should be generally lighter 
than W- or Z-boson. In this case, gravitino mainly decays into photon and neutrino. The strength 
of the 7 + z/ decay that occurs through the diagram depends on the value of the photino-neutrino 
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mixing. In our model this mixing occurs via the oneloop diagram in Fig. 4 and the corresponding 
mixing parameter Ux^-v can be estimated to be 

_ < > v^ksinpfi'^ehe 

and is of order 10"^ - 10"^ for jj. ~ MgusY ~ ~< >~ 100 GeV and K ~ 10"^ The 
decay rate is then given by 



r(G'^^ + 7)«i^^, (20) 

^ ^' 327r M2 ' ^ ^ 



where iC/i^-yl is the percentage of photino in neutrino mass eigenstate in the presence of R-parity 
violation v ^ Ui,^^ H as in Fig. 4. 



<V^> 




®H 

U 



FIG. 4: 1-loop diagram for photino and neutrino mixing. 



The transverse and longitudinal components of a decaying gravitino give comparable couplings 
to final state photon and neutrino. As a rough estimate, taking 7723/2 ~ 2 GeV and \Ui,;y\ ~ 
10-'5- 10-7 (as given in our model), one gets 1(6" ^ v+iY^ ^ lO^^-lO^^ GeV"^ ^ lO^^-lO^^ 
sec. Thus we see that even though the gravitino is unstable via R-parity violating interactions, it 
can become a viable dark matter of the Universe. The gravitino also has three-body decay modes 

with r3_^^^y = I ig2n'iMpiMt — ) ^^^^ ^^^^ small compared to the two body decay rate. 

The gravitino decay also produces energetic extra-galatic diffuse 7-rays in the universe. From 
the above decay rate, one can derive the photon energy-flux per unit solid angle assuming gravitino 
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has homogeneous distribution throughout the universe 121] 



dE 



Pef^3/2r(G'^ I/ + 7) / 2E 



87rr3/2-ffofi 



1/2 
M 



2E ^ 



5/2 



'1/2 



"^3/2 



1 + 



2E 



-1/2 



(21) 



where the last step follows if G ^ z/ + 7 is the dominant decay channel of gravitino. A ^^-function 
9(1 ) is implicit because the photon frequency is cut off at 7713/2/2, where the flux is 

V "^3/2/ 

peaked. Experimentally, EGRET observes the extragalatic cosmic gamma-ray flux has an excess 
over the power law spectrum, which can be as large as 2.23 x 10^^ (cm^ str s)^^ GeV for photon 
energy between 2 and 20 GeV ||24|]. Attributing part of this excess to the decay of gravitino dark 
matter in the universe, one has an upper bound on the mass of gravitino for a given neutrino- 
photino mixing \U„y\. Taking the mixing parameter around 10^^ we find that 7713/2 ^ 2 GeV for it 
to be consistent with cosmology. 



VI. NLSP AND VERTICED DISPLACED AT LHC 

The above R-parity breaking scenario can be tested at LHC. Here we discuss the possible 
signatures following from various decays of various possible NLSP's. 

A. Neutralino as NLSP 



R-parity conserving decays The next-to-lightest superparticle can be the neutralino, i.e. a linear 
combination of neutral gauginos and higgsinos. It has both RFC and PRV decay channels. The 
R-parity conserving channels are^* — > G + 'j,G + Z'^,G + h^, with gravitino dominantly in the 
longitudinal component. These decay rates are calculated in the same way to gravitino decaying 
into neutrino and photon yielding 

|f/5^0^pp777|o 



r(x°- 



r(x° - Gh' 



167rF2 
167rF2 



(22) 
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If/^fOjp is the percentage of the i-th species (photino, zino, higgsino) in the neutralino NLSP and 
"I is the coupling constant of goldstino with matter fields. As will be seen, such decays are much 
slower than R-parity violating modes. 

R-parity violating two-body decays The light neutral gaugino (wino or bino) does not directly 
couple to right-handed sneutrino (standard-model charge-free), so the higgsino component in 
NLSP controls its R-parity violating decays. If the neutralino NLSP is heavier than the W or 
Z-boson, it can decay into + or + u. The corresponding Feynman-diagrams are shown 
in Fig. [51 and the decay rates are on the order of 




FIG. 5: A typical two-body (left) and three-body (right) PRV decay diagram of neutralino NLSP. 

To make a rough estimate, we take \U^Oy\ = 10~^ and m^o = 100 — 200 GeV. We then find 
that r2_body ~ 10"^^ — 10"^^ GeV which corresponds to a neutralino lifetime of ~ 10"^^ — 10"^^ 
sec. This will lead to a vertex displacement in the detector of about 0.1 — 1 mm. 

This decay rate is fast enough not to ruin the success of BBN, and slow enough to produce 
collider signatures such as vertex displacement at LHC. The displacement of secondary vertex 
where NLSP decays is around cr ^ 1 mm, which is observable within the detector. 

R-parity violating three-body decay If the neutralino is lighter than + ox + u mass, 
the two-body RPV decay channels are forbidden. It can only decay into a three-body final state 
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60 80 100 120 140 160 180 200 
(GeV) 



FIG. 6: Decay rates as a function of neutralino mass. Red curve is for — > W^l^ , blue is for ^ Z^v, 
and green is for ^ —>■ ULbihR, where we choose \U^y\ k, 10^^ and \U^^\ ~ 0(1). We also choose 
tan (3 = 2) and rriy = 300 GeV. 



through a virtual W- or Z-boson or sparticles, which couples to two SM particles. The dominant 
three-body decay mode is going to h + h + v [25], as shown in Fig.[5l The three-body decay rates 
have been calculated in Ref. UM- In Fig. 6, we plot the two- and three-body decay rates as a 
function of the neutralino NLSP mass. We see that the three-body decay has a longer lift time and 
can hardly produce an observable vertex displacement inside the detector near the vertex. 



B. Stau as NLSP 

The stau can also serve as the NLSP in the model. Generally, there is a small mixing between 
left- and right-handed staus. The lighter stau ri is a linear combination of them 



(24) 



where |ap + = 1. From the discussions in Section IV, we know that at the R-parity violating 
vertices, the right-handed sleptons mainly couple to lepton final states, while left-handed sleptons 
can decay either leptonically or hadronically. Typical diagrams for stau decay are shown in Fig. |71 
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The decay rates are 



T{tl ^ ih) 
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(25) 



where the mixing due to RPV mass-insertion can be estimated to be \Ur^H- 1 , \Ui^ho i 



10- 




<v^> 




FIG. 7: Feynman diagrams for stau decay. 
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FIG. 8: Stau NLSP decay rate as a function of its mass. Blue curve is for |ap = 0.05 and the red one is for 
|ap = 0.3. There is a kink at = rrit where a significant hadronic final-state decay channel opens. 



In Fig. 8, we plot the stau decay rate as a function of its mass. There is a kink at m:^^ = 
where a significant hadronic final state decay channel opens. For stau mass less than 300 GeV, the 
decay rate is of the order of lO^^'^ GeV, which corresponds to a vertex displacement 1 mm. 

Actually, the NLSP can also be left-handed sneutrino. The decay diagrams are similar to those 
of stau NLSP, and one the sneutrino has a similar decay rate to sleptons. 

18 



VII. ADDITIONAL COMMENTS 



In this section, we discuss some other aspects of the model. 

Strong and SUSY CP It has been pointed out ll29ll that the constraint of left-right symmetry 
restricts the mass matrices and phases in the model in such a way that it provides a solution to both 
the strong CP and SUSY CP problem. Key to solving the strong CP problem is the hermiticity of 
the quark mass matrices in the model. Once R-parity is broken, it is not clear that this property 
will hold. However we find that if the right handed sneutrino vev is aligned along the electron 
flavor, the induced phases in the Det are of order (^^^7^)^ ~ 10^^° which is below the bound 
on the 9 provided by electric dipole moment of the neutron. SUSY CP problems are not affected 
by this phase. 

Decays of the doubly-charged Higgsino One of the distinguishing features of the B-L=2 
triplet Higgses is the presence of doubly charged Higgs fields. These are accompanied by their 
fermionic superpartners (A++). These particles can have masses in the sub-TeV to TeV range 
and may therefore be accessible to LHC. In the presence of R-parity breaking, these fermions will 
decay to A++ — > rHb (as shown in Fig. [9l). There are no standard model background for such 
decay modes. 




FIG. 9: A Feynman diagram for doubly-charged Higgsino decay. 



R-parity breaking processes in the early universe We now make a few more comments on 
the impact of R-parity violating in our model in the early universe. Note that due to the smallness 
of the strength of the eiLHu term in the superpotential i.e. Cj ~ 10~^ GeV, the effective A {LLe^) 
and A' {QLdf) interactions are of order or less than equal to 10^^. It was noted in ref.[28], for 
such small strengths, the lepton number violating interactions involving sparticles are out of equi- 
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librium. As a result one could perhaps contemplate generating lepton asymmetry by the decays of 
the NLSP particle e.g. neutralino slightly above the electroweak scale and have them converted to 
baryons via the sphaleron effects. This question is under consideration. 



VIII. CONCLUSION 



In conclusion, we have discussed the implications of the minimal renormalizable supersym- 
metric left-right seesaw model as an extension of MSSM to include neutrino masses and show 
that the model leads to dynamical R-parity violation by its ground state in order to break parity 
symmetry. We have analyzed a particular realization of the model where spontaneous R-parity 
breaking occurs along the u^e direction. We first show that small neutrino masses in this model 
can be understood only by the usual low scale seesaw condition that the Dirac Yukawa couplings 
are tuned to the value of < 10^^ and without any further tuning. We then show that if gravitino 
is the LSP, then its decays are automatically suppressed by the same condition that guarantees the 
smallness of neutrino masses, making the gravitino long lived enough to be an unstable dark matter 
of the Universe. We also point out that the NLSP decays in this model may lead to displaced ver- 
tices which can provide a clear LHC signal. These models have many other collider implications 



such as doubly charged Higgs and Higgsino fields BOh that have been discussed extensively in the 



literature as well as low energy lepton number violating signals such as muonium-anti-muonium 
oscillations Ibi]]. 
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Appendix 

In this appendix, we discuss the minimization of the Higgs potential to obtain an upper bound 
on the right-handed scale vr. 

First we will show that in the SUSYLR model, if R-parity is preserved by the vacuum, parity 
symmetry cannot be broken i.e. vr = vr = Q d& shown in [8J. A more precise statement perhaps 
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is that, if we express the potential, V as functions of the vevs of the neutral Higgs fields ( vr, vr 
and X and look for a minimum of V along the direction x = (i.e. R-parity conserving ), the 
minimum occurs at vr = vr = 0. We then show that once we include R-parity breaking effect by 
the vacuum i.e. (z/^) = x ^ i.e. along the direction x ^ 0, there appear global minima that break 
parity i.e. vr, vr ^ and also that it occurs only below a certain value for vr and vr i.e. there 
is an upper limit on the parity breaking scale. This proves our assertion that R-parity breaking in 
this theory is a dynamical phenomenon. 

To show this let us start with the potential in Eq. (2) which consists of the field vevs vr and vr 
and X = {u^) and look for its minimum: 



V = [MIvI + MIvI~2\B\vrVr] 

+ [fx^ - {2AvR + 2f^i^VR -ml- Afvl)x^] + 



(26) 



We have set the $ and A, A vevs to zero. Note that we have kept the < z/^ > in the potential. For 
simplicity, we have set / = /i. 

The first point to note is that to ensure a lower bound on the potential, we must satisfy the 
conditions on the parameters: 

Ml + M\>2B , 

MaM^ > B . (27) 

The first constraint comes from looking at the direction vr = vr and x = demanding that 
the potential is bounded from below. The second comes from looking along the QED breaking 
vacuum so that the D-terms vanish and setting x = and again demanding positivity. Once these 
two conditions are imposed, for x = 0, the minimum of the potential corresponds to vr = vr = 
and hence no parity violation. 

It is worth pointing out that the form of the potential for x = is same as in the case of MSSM, 
where of course we know that symmetry breaking occurs. The difference in the case of SUSYLR 
is the observation[|8D that there exist QED breaking directions along which for arbitrary vr and vr, 
the D-term vanishes so that one has the second condition in Eq. (|27l) . In the case of MSSM, the 
second condition does not exist and in fact to break the gauge symmetry in SUSYLR case, one 
needs the opposite of the second condition i.e. M^M^ < B. 
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In order to show that along the nonzero x directions, one can indeed have a minimum that 
breaks the gauge symmetry, it is convenient to rewrite the potential in Eq. (|26|) as follows: 



V 



C2 , 9\,2 



{f + '-)[x'-CY + D{vn,vn) 



(28) 



where 



D{vr,Vr) 



c 



Mlvl + Mlvl -2\B\vnVR + 



C2 , ~9' 



[2Avr + 2f^l^VR -ml- Apvl + \g\v 



R 



(29) 



2(p + i) 

where we define = {gl + g"^). Advantage of rewriting this way is that we can now minimize 



with respect to x very easily and get 



X 



C. 



(30) 



Since x is a real number, the above equation implies that C > 0. For vr and vr close to each 
other, the C > condition turns into an upper limit on the vr scale of 



{A + //iA) - ^[A + fii^Y-APml 



<VR< 



{A + ffi^) + V(^ + //iA)2-4/2 



. (31) 



4/2 - 4/2 

This is clear from the expression for C since for large vr, the —pvl term in C dominates mak- 
ing C < and hence driving the vev of x to zero in which case, the minimum corresponds to 
Vr = Vr = as noted and hence no parity violation. This is in accord with the observation 
of ref.[|9(] that the tree level potential for the minimal SUSYLR model requires R-parity violation 
if parity has to break, as required to get the standard model and the parity breaking scale has an 
upper limit in the TeV range, making the theory experimentally testable at LHC. 

We have done a numerical analysis of the potential and find that the full potential indeed has a 
negative minimum value only when the x ^ and vr ^ and vr ^ giving the desired parity 
violating graound state. We also find from this numerical analysis that once we set C < or 
equivalently x = 0, the ground state corresponds to vr = vr = 0. In the table below, we 
give some numerical examples of solutions for the desired vacua for specific TeV scale parameters 
in the potential for the case when x ^ 0. We check that the potential has a minimum for each 
case with a value of the potential at the minimum which is negative so that it is indeed a global 
minimum. 
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/ 


A (GeV) 


^ (GeV) 


mo (GeV) 


VR (GeV) 


VR (GeV) 


(GeV) 


0.1 


100 


200 


500 


2363 


1677 


2286 


0.1 


100 


500 


600 


2897 


2166 


2658 


0.2 


500 


200 


500 


3561 


2304 


4020 


0.3 


500 


500 


500 


1865 


1736 


1815 



TABLE I: The other parameters are chosen to be A/a = 1 TeV, = 1.1 TeV and B = 1.099 TeV^ 



Incidentally, along the x = direction in the potential, there is also no electrweak symmetry 
breaking by similar arguments as above [|8D i.e. (0) = 0; however along the above x 7^ direction, 
one immediately gets (</>) 7^ also along with breaking of SU{2)r x U{1)b-l symmetry. 
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